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Abstract 

A probability method is provided to prove three classes of com- 
binatorial identities. The method is extremely simple, only one step 
after the proper probability setup. 

1 Introduction 

The probability method has been a popular and powerful tool to solve dis- 
crete problems in the last half century since Paul Erdos. It creates a new 
branch of mathematical research pQ. One of the beauties of this method is 
that there is being creative and open-minded involved, especially when the 
discrete structure and the probabilistic construction look totally irrelevant 
at first sight. This paper gives a probability method to prove some combi- 
natorial identities. 

The first result is a class of combinatorial identities involving symmetric 
functions, 

n 

^2(xj) m Y[(xj - Xi)- 1 = h m _ n+ i(x), (1.1) 

J=l i¥=3 

where m G Z + U {0} and h m _ n+ i(x) is the (m — n + l)-th homogeneous 
symmetric function of X\, x n . 
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The combinatorial methods were provided by Louck and Biedenharn [9], 
Strehl and Wilf [32] and Chen and Louck [3]. Equation (jl.ip is a generaliza- 
tion of Good's identity [7], 
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The summation on the left hand side of (11. ip is important in the explicit 
matrix elements in the unitary groups, as well as multivariable hypergeomet- 
ric series well poised in SU(n). More combinatorial and algebraic applica- 
tions of (II .ip and (11.21) can be found in [5] and [PJ. 

The second result contains multinomial convolutions. 
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where m > 1, > 0, > and A = J^ILi a i- The {£.,■} are defined recur- 
sively in Section 3. 

A special case is about the convolution of two sequences. 



n „ oo n 

j=0 i=l 



The combinatorial proof of this type of identities was discussed in [11] and [8] . 
The third result gives an identity involving multinomial convolution too. 
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where A = YH=i a i an d 
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A similar equation of the same type involving Bernoulli numbers was 
proved by Dilcher [I]. 

A special case provides another method to prove an identity about the 
Stirling numbers of the second kind. 



This probability method takes advantage of the nice structure of some 
probability density functions. In the next three sections, the combinato- 
rial identities are proved by the summation of exponential random variables, 
Gamma random variables and uniform random variables. The last section 
discusses some inherent connections between the probability density func- 
tions and the combinatorial structures. 

2 Exponential Distribution 

The probability density function of an exponential distribution with param- 
eter A > is 



In this paper, an important component of the proof is the density function 
of the summation of independent random variables of some specific distribu- 
tion. The probability density function of the sum of independent exponential 
random variables with distinct parameters is as follows, which can be found 
in Feller [5], #12, page 40 (with no proof). 





Lemma 2.1. Let {Xi, X n } (n > 2) be a sequence of independent ex- 
ponentially distributed random variables with distinct parameters Ai,...,A n 
(Aj > 0), respectively. The density function of their summation YTj=i^j ^ s 



n n 
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fn(x)=[H\ J ]j2^ , ^>0- (2.1) 



1=1 
l^k 



Proof: When n = 2, by convolution 
f 2 (x) = AiA 2 



A2 — Ai 

Suppose it is true for n — 1. Then from convolution, we can get 



/n(*)=[n^]£ 

3=1 k=l 



1 II(A*-A fc ) 



z=i 



In order to finish the proof, we need to show that 



n—l 1 

-Y — - 



n-1 



=1 II( A '-^) H(A*-A„) 



Or equivalently, 



=1 

l+k 



- 1 n< A ' - a '» 

1=1 

Multiplying the common denominator on both sides, this is equivalent to 



£(-l) fc II (Ai-A i ) = (2.2) 



k=l l<j<l<n 
l+k 
3^k 
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Note that the product is a Vandermonde determinant: 
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From this, it is not hard to see that (12. 2p is the determinant of 

/ 1 i x 1 ... Ar 2 \ 



1 1 A fc 
V 1 1 A n 



Ar 2 



Ar 2 j 



which is obviously zero. 



Theorem 2.2. For any distinct positive real values Ai,...,A n (n > 2) and 

nonnegative integer m, 

Xi ^ 1 



\m W 



k=1 K ¥k A/ - Xk 
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Proof Let {Xj}, 1 < j < n, be independent exponential random variables 
with parameter Xj, respectively. The m-th moment of is 
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On the other hand, by integrating the density function (12. ip . we have 



n n . 

«[e>r]=££n 
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Equation (11.11) is an easy consequence of Theorem 2.2. 
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3 Gamma Distribution 

A Gamma random variable X has probability density function 



e 



-x//3 



f( x ) = x a ~ x ~, — , x>0, 
J v ; T(a)ft a ~ 

where a > and ft > are two parameters. 

The density function of the summation of independent Gamma random 
variables is given in [TO] . 



Lemma 3.1. Let Xx, ■ • • , X n be independent Gamma random variables. The 
probability density function of Xi is 



with aci > and /3j > 0. Without loss of generality, we may assume that 
ftx = min {fti}. 

l<i<n 

Let S n = Yli=i Xi ■ Then the density function of S n is 

9n ( x ) = rr (ft.) Ta./+m — e - — x > o 

ii w U nA +J )ftt +j 

where A = J^ILi a « anc ^ ' s are defined by the recursion 

j+i 



1 

= — T Yl l ^l 5 3+l-h j = 0, 1, 2, 

So — 1 and 

n 



i=i 



Using two methods to compute the expectation of S n , we can get an 
identity involving the convolution of the two sequences {c^} and {fti}. 
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Theorem 3.2. 



n „ oo n 

i=l ^ % j=0 i=l 
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with A and Sj as defined in Lemma \3. 11 

From the computation of moments, a general identity is derived. 
Theorem 3.3. If m is a positive integer, 



4 Uniform Distribution 

A uniformly distributed random variable X on the interval [a, b] has proba- 
bility density function 



Two types of summations are considered. One is the summation of I.I.D. 
(independent and identically distributed) uniform random variables. The 
other one is the summation of independent random variables uniformly dis- 
tributed on different intervals. Before the discussion, the following notation 
is needed. 



Let Xi, X n be uniformly distributed I.I.D. random variables. Without 
loss of generality, we may assume that Xi is uniform on [0,a]. (If X^ is 
uniformly distributed on [b, c], then let Yj = Xi — b. So Y^s are I.I.D. random 
variables uniformly distributed on [0, c — b] and J27=i = SILi ^ + 
And it is sufficient to consider Y^7=i The following lemma, which can be 
found in [5] page 27, gives the density function of the summation. 





x + \x 



, for iGR. 
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Lemma 4.1. For n> 1, let S n = Y^i=i -^-i- The probability density function 
of S n is 

M^s^nyi D-d 1 (")(*- »)r'. ("J 

/on e [0, na] and = otherwise. 

Using two methods to compute the m-th moment of we can get 
Lemma 4.2. For m > 1, 

]_ , n m+n /^X Al — 1\ _ jm+j+l-i . 

2 — 1 J — U 

= V f m ^ - (42) 

n + .^t= m vi,---,W(^i + i)---(^ + i) 

The right side of the above identity can be simplified. 
Lemma 4.3. 

x - / m \ 1 S(m + n,n) 

U • • • , < J (ii+i)-(in+i) = ( m r) ' 

where S(m + n, n) £/ie Stirling number of the second kind. 
Proof: Let 



1m, n ^ ^ 



m 



>W (ii + 1) • • • (i n + 1)' 



iiH Mn=m 

Note that 7 mj „ is the coefficient of t m+n /m\ in (e* — l) n . Hence, 



,»(e*-i) = y L.( e *-ir 

^ n! 

n>0 

„.n 4-m+n 



„! 



n! 14 — ' ml 

n>0 m>0 



"7tn,n 
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Besides, 



n>0 k>n 



X> B [$>(ro + 



j-m+n 

n, n) 



(m + n)\ 

n>0 m>0 v ' 



Comparing the coefficients, there is 

— r 7m,n = y^S(m + n,n) 

n.< z — ' ml z — ' 



So 



n! ^— ' m! ' ^— ' (m + n)! 

m>0 m>0 v y 



D(m + n, n)- 



n! (m + n)! 



Observe the density function (14.11) . The m-th moment of S n can be 
represented in another way, instead of the two sides in (14. 2p . 

Lemma 4.4. 



Proof: First, 

rn 

(x-i) n ~ l x m dx = / (x - i) n - 1 x m dx 



■m+n 



m!(n-l)!yV ]\j/ m + n 

(m + n)\ +^ \m — j 

m\{n-l)\ S^ ( _ ir+j ( m + n 

(m + n)! \m-j 



[n - i) n+j n m - j 



(4.3) 
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Note that the summation in (14.31) is part of the expansion of i 



m+n 



m+n 



.• ( m + n 



3=0 



i m+n = "' . " \{n - i)'n"-''--> 



J 



But by symmetry, the other part will vanish after the summation of i. 

n— 1 



i=0 



3=0 



m + n 



3 



(n - i) j n m+n - j 







Hence, the uniform density function gives another method to prove the 
identity of Stirling function of the second kind. 

Theorem 4.5. 



S(m + n, n) 



' i=0 ^ ' 



If Xi , . . . , X n are independent random variables and Xi is uniformly dis- 
tributed on [0, dj], the density function of their summation is given by [12J. 

Lemma 4.6. For n > 2, let S n = Y^=i The probability density function 
of S n is 



h n (x) 



[n 



l)lJJ ai " i=1 l<ii<J2<-<j!<n 



i=i 



n-l 



Theorem 4.7. For m > 1, 

1 r _4 m +™ 



m + n 



Bin) 



E 



n H Hn=m 



m 



k, ■ ■ ■ , In] k + 1 «n + 1 



i=l 



where A = Yli=i a « anc ^ 



B(n) = B-D' E [E 

i=l l<ji<j'2<---<ji<n k=0 



n-l 



n — 1 



E 

i=l 



ri— 1— k 



k+m+l 



1=1 



k + m + l 
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5 Relations and Future Development 



It is easy to start with the summation of some independent random variables 
to get some identities, if the density function is provided. But given an 
identity, how to construct a probability distribution for the proof is another 
question. Hence, it is important to find out the inherent relations between the 
probability density functions and the combinatorial structures. One possible 
connection is their "generating functions" . All of the identities in this paper 
were studied by the generating function method before. The generating 
function of a probability distribution is its characteristic function (ch.f.). 
The ch.f. of a random variable is defined as 

<p x (t)=He ltX ). 

For example, the ch.f. of an exponential distribution with parameter A is 

m 1 



i it 

1 A 



And Lemma 2.1 can also be proved by the ch.f. 4> n (t) of the summation 

n 



3=1 



n 

, n (t) = l[Ee UX ' 



3=1 



n 



3=1 \ 

Note that this characteristic function has the same form as the generating 
function of the homogeneous symmetric functions h r (x) in variables x±, X2, ■■■ 

G(t)=J2hX = H(l-x ] t)- 1 . 

r>0 j>l 

Besides, the density function (12.11) has similar product inside the summa- 
tion to the form of Lagrange interpolation polynomial 

n 

3=1 k^j X 3 Xk 
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In fact, Good's identity (ll.2p and its generalization fll.il) can also be 
proved by Lagrange interpolation method (as in [3] and [T3]). 

If we compute the expectation E(exp(t ^2 Xk)) by integrating the integral 
transform of (JUD and by the fact that E(exp(i^X fc )) = ]J E(exp(tX k )), 
we get another identity 

n . n 

T^ — Tlx^TrR— t<min{X u ...,X n} . (5.2) 

k=l A fe i-Lk k=l 

Observe that (I5.2p can deduce the interpolation proof of Theorem 1 in [3] 
and [13]. 



We may also check the moment generating function (m.g.f.) of a distri- 
bution, if it exists. The m.g.f. of a random variable X is defined as 

M x (t)=E(e tx ). 

The m.g.f. also has the property that the m.g.f. of the summation of inde- 
pendent random variables is the product of their m.g.f. 's. And the moments 
can be got from the m.g.f.. 

E(X m ) = M ( x\0). 
For example, the m.g.f. of the uniform distribution on the interval [0, 1] 

is 

e* - 1 
M(t) = C - T ±. 

Then the m.g.f. of S n , the summation of n independent uniform random 
variables on [0, 1], is 

(e f - l) n 

Note that the exponential generating function of the Stirling number of 
the second kind has similar form. 

Fn(t) = S{m + n)- - 

w ^ v ' (m + n)\ 

(e* - l) n 
n\ 
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Look at the exponential generating function of ( m ^ n )E(5 1 ™). 



z — / \ n I (m + ny n\ z — ' m 

m>0 v 7 V 1 7 m > 

T (e* - 1)" (e* - 1) 



n\ t n n\ 
This proofs Lemma 4.3. 
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